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Abstract 



We describe hot, optically-thin solutions for one-temperature accretion disks around black 
holes. We include cooling by synchrotron, bremsstrahlung, and Comptonization. Our so- 
lutions are thermally and viscously stable, with gas temperatures on the order of T ~ 
10 9 - 10 10 7 K. The thermal stability is a direct result of the inclusion of synchrotron cooling. 
The new solution branch is related to the advection-dominated solution for a two-temperature 
, gas described by Narayan & Yi (1995b). It is present only for mass accretion rates less than 

' some critical Merit which depends on the radius R and viscosity parameter a. The solutions 

are advection-dominated for extremely low values of M. However, for a range of interme- 
diate accretion rates, the new solutions are both hot (T ~ 10 10 K) and cooling-dominated. 
Because of this new feature, one-temperature solutions are significantly more luminous than 
the corresponding two temperature solutions. 
\& , The radial profile of the new solutions is unusual. The inner parts of the flow are cooling- 

dominated and have a disk-like geometry, while the outer parts are fully advection-dominated 
^ c"| and nearly quasi-spherical. 

P 1 Introduction 

£3 ■ 

The best known model of an accretion flow is the standard thin accretion disk developed by 
^ ! Shakura & Sunyaev (1973), Novikov & Thorne (1973), and Lynden-Bell & Pringle (1974) (see 

ITji \ Frank, King Sz Raine (1992) for a detailed discussion). The basic assumption of the model is 

that the accreting gas is cool, compared to the local virial temperature, so that the flow acquires 



a thin disk configuration. The disk parameters are calculated assuming an equilibrium between 
the viscous energy generation inside the disk and radiative cooling from the surface; the latter is 
computed under the assumption that the disk is optically thick in the vertical direction. 

The thin disk model has been widely and successfully used for modeling various low energy 
systems such as accreting white dwarfs and pre-main-sequence stars (Frank et al. 1992). However, 
the model has had less success explaining the characteristics of relativistic objects such as accret- 
ing black holes. There are essentially two problems: 1) In a system containing an accreting black 
hole, the accretion disk must extend inwards to the Schwarzschild radius, Rschw However, for 
reasonable mass accretion rates and the standard a viscosity prescription, the Lightman-Eardley 
instability (Lightman & Eardley 1974) causes a breakdown of the thin disk configuration at radii 
R ~ Rschw Therefore, there is some doubt that a thin accretion disk can exist at all close to 
a black hole. 2) Even if a thin disk is viable, there is a problem with the spectrum. Since the 
gas in a thin disk model is optically thick, we would expect a roughly blackbody-like spectrum 
at the high frequency end, with a cutoff at several keV. However, extensive observations of black 
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hole systems by GINGA, GRANAT, and GRO (Tanaka 1989, Grebenev et al. 1993, Maisack et 
al. 1993, Harmon et al. 1994, Johnson et al. 1994, Kinzer et al. 1994, Gilfanov et al. 1995) have 
made it clear that the spectra of most black hole accretors have a hard power-law components 
extending to few x 100 keV. This component has to be radiated by an optically thin plasma with 
temperature T ^ 10 9 K. The required temperature is at least an order of magnitude higher than 
that predicted by the standard thin disk model. 

Shapiro, Lightman, & Eardley (1976, hereafter SLE) discovered a new class of solutions at sub- 
Eddington accretion rates, where the accretion flow is optically thin and quite hot (T e ~ 10 8 — 10 9 
K). These authors also introduced the important idea of a two-temperature plasma, in which ions 
are much hotter than electrons. Since the SLE solution has exactly the necessary characteristics 
(optically thin hot gas) to explain the spectra of accreting black holes, it has been widely studied 
and applied in models of X-ray binaries and active galactic nuclei (e.g. Kusunose & Takahara 
1985, 1989; White & Lightman 1989; Wandel & Liang 1991; Melia k Misra 1993; Luo & Liang 
1994). Unfortunately, the SLE solution is thermally unstable (Pringle, Rees, & Pacholczyk 1973, 
Piran 1978). The efficiency of bremsstrahlung cooling in the optically thin gas is proportional 
to the particle density which decreases with increasing temperature. Therefore, as the gas is 
perturbed to a higher temperature, the rate of cooling per unit mass decreases and the gas heats 
up even further in a runaway process. The thermal instability of the SLE solution makes it 
unlikely that real flows can take up this configuration. 

Recently, a new class of two-temperature advection-dominated solutions has been discovered 
(Narayan & Yi 1994, 1995a, 1995b; Abramowicz et al. 1995; Chen 1995; Chen et al. 1995). 
These advection-dominated solutions are optically thin, are hotter even than the SLE solution, 
and are viscously and thermally stable to large-wavelength perturbations (Abramowicz et al. 
1995, Narayan & Yi 1995b). Although Kato, Abramowicz & Chen (1995) have discovered an 
instability in these flows at short wavelengths, they show that the mode amplitude does not grow 
significantly and so the viability of the solutions is not in doubt. The advection-dominated model 
has been fairly successful in explaining a number of low-luminosity systems (e.g. a model for 
Sagittarius A* by Narayan, Yi, & Mahadevan [1995], a model for the soft X-ray transient A0620- 
00 by Narayan, McClintock & Yi [1996], and a model for the central source in NGC 4258 by 
Lasota et al. [1996]) and is perhaps also a reasonable model for more luminous systems (Narayan 
1996; Narayan, Yi, & Mahadevan 1996). 

The most detailed advection-dominated models considered so far are based on a two-temperature 
plasma which cools via bremsstrahlung, synchrotron and Comptonization processes. However, 
there have been discussions in the literature questioning whether a two temperature plasma can 
occur at all in astrophysical accretion flows (Phinney 1981, Rees et al. 1982). The discussion 
has been further fueled by the work of Begelman & Chiueh (1988), who identified a particular 
mechanism involving collective plasma waves, which may under certain conditions pump en- 
ergy directly into the electrons and bring the ions and electrons into thermal equilibrium. An 
important question therefore is the following: are there stable hot solutions, analogous to the 
two-temperature advection-dominated solution, if the plasma is well-coupled and has a single 
temperature for the ions and electrons? Abramowicz et al. (1995) and Chen (1995) did consider 
single-temperature plasmas in their studies of advection-dominated flows, but they included only 
bremsstrahlung cooling. As yet, one-temperature models with a more detailed cooling, including 
synchrotron radiation and Comptonization of synchrotron and bremsstrahlung photons, has not 
been considered, even though at T ~ 10 10 K with an equipartition magnetic field these processes 
completely dominate over pure bremsstrahlung cooling. This is the study we present here. 

In §2 we derive the equations that describe an accreting one-temperature plasma. We then 
present the main results of the paper in §3. We find that, as in the two-temperature case (Narayan 
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& Yi 1995b, hereafter NY), the equations allow three branches of solutions: the standard thin 
Shakura & Sunyaev disk, a hot optically-thin, thermally unstable solution, which is equivalent to 
the SLE branch for a two-temperature plasma, and finally a new hotter and thermally stable disk 
solution. We describe the properties of the new solution and compare it to the corresponding two- 
temperature solution branch. Among other things, we find that the one-temperature solution is 
hotter than the equivalent two-temperature one and is not advection-dominated for a wide range 
of accretion rates. This means that for a given M the one-temperature disk is considerably more 
luminous. We conclude in §4 with a summary and discussion. 



2 Basic Equations 

2.1 One-Temperature Accretion Flow 

Narayan & Yi (1994, 1995b) derived a set of local equations describing a vertically averaged, 
axisymmetric, two-temperature advection-dominated flow. Here we slightly modify these equa- 
tions by setting the electron and ion temperatures equal to each other. We also include radiation 
pressure, which NY neglected. 

The vertically averaged pressure p and density p of the accreting gas are given by: 

where c s is the local isothermal sound speed, M is the mass accretion rate, R is the radius, H is 
the vertical scale height of the disk, and v is the radial velocity. We write the total pressure as the 
sum of the gas, radiation and magnetic pressures, p = p g +p r +Pm, and we define three quantities 
that describe their relative magnitudes: (3 = p g /p determines the fraction of the total pressure 
due to gas pressure; (3 m = (p r + p g )/p = 1 — p m /p determines the importance of magnetic field 
pressure; finally % is defined such that = xflm so that % represents the ratio of gas pressure 
to the sum of gas and radiation pressure, x = Pg/(Pg + Pr)- For x ~ 1 radiation pressure is 
unimportant, and for \ ^ 1 radiation pressure dominates over gas pressure. In our calculations 
we fix the value of f3 m , typically at 0.5 corresponding to magnetic pressure equal to 50% of the 
total pressure, and then solve for x and /3. 

With the above definitions, the scaled relations derived by NY take the form: 



v 



-2.12 x 10 10 ^ cms" 1 . 



IT 



c 2 = 4.5xl0 20 ^cm 2 s- 2 , 



r 



p = 3.79 x 10~ 5 m / gem" 3 , 



B = 6.55xW ^ ^G, 



ac\m V r<- 



f = 1 p/3- 7 \ ^ 32 - 24/3-3/3 5 



7-1 ) ' ' 24-21/? ' 

q + = 1.84 x 10 21 em y^ erg cm -3 a" 1 . (2) 

Here B is the magnetic field strength and the quantity / determines the degree to which the flow 
is advection-dominated. / ranges from to 1; the limit of / = 1 corresponds to the extreme case 
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when all the energy is stored in the gas and advected, while / — * corresponds to a standard 
cooling-dominated disk where most of the generated energy is radiated locally. The quantity q + 
is the viscous dissipation of energy per unit volume, and the factors e',ci, and c 3 are given by 



e , 1 I 18a 2 -x 2 c \ _ 9c 3 a 2 



2 \ 2x 



ci = -c 3 



3 

Tr 



c 3 = 3.12 x 1(T 13 — . (3) 

In the above equations masses are scaled in units of the solar mass 

M = mM Q ; 

accretion rates in Eddington units 

M = mM Edd , Mem = = M = 1-39 x 10 18 m gs" 1 , 

where K es = 0.4 cm 2 g _1 and we assume the standard efficiency factor r\ e ff = 0.1 (e.g. Frank et 
al. 1992); and radii in Schwarzschild units 

2GM . 
R = rRschw, Rschw = — o— = 2.95 x 10 d m cm. 

<r 

If we choose the values of r, m, rh, a and (3 m , and assume specific values for (3 and T, 
equations ([[]), @ and (|3|) yield all the other disk parameters. 



2.2 Electron-Positron Pair Equilibrium 

Some of the hot solutions we calculate in this paper have temperatures of up to 10 10 — 10 11 K (see 
Figure §[d]), and it therefore becomes necessary to include in our model the effects of relativistic 
pair production and annihilation. Since a complete solution of the pair balance problem is beyond 
the scope of our calculation, we make some simplifying assumptions. 

Svensson (1982, 1984) has explored the properties of pair equilibria in a uniform relativistic 
plasma where photons are generated via pair annihilation and bremsstrahlung. He finds that for 
temperatures below a certain critical value T c there are two equilibrium branches available that 
are characterized by the ratio z of the pair number density, n+, to the number density of protons, 
n p = n e — n + : (1) an optically thin low-z branch {z <C 1) where pair production is dominated 
by particle-particle processes, and (2) a high-z branch [z > 1) where the Thomson scattering 
optical depth is of order unity and photon-photon pair production dominates. We confine our 
calculations to the low-z branch where only particle-particle processes are important. Svensson 
has shown that the electron-electron pair production rate is always larger than the electron- 
proton rate by a factor of ~ 10. This allows us to ignore the contribution from electron-proton 
collisions. 

With the above simplifications, we can solve for the equilibrium pair density, z, analytically at 
any given temperature. For a thermal distribution of electrons and positrons, the pair annihilation 
rate per unit volume is (Svensson 1982): 

(n+)ann = ncr 2 n e n + g(9) , (4) 
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where 9 = kBT/m e c 2 is the dimensionless electron temperature, r e = e 2 /mc 2 , and g(9) is given 
by: 



9(0) 



1 + 



(5) 



In (1.120 + 1.3) 

For the electron-electron pair production rate we adopt the expression used by White & Lightman 
(1989): 

,. , 2 2 J 2 x !0~ 4 3/2 exp (-2/9) (1 + 0.0150), if < 1, 

(n+j ee - cr e n e j (m/277r) ft 2 ^3 (1 + O .O58/0)-\ if 9 » 1, W 

where aj is the fine structure constant. Local equilibrium requires that (n + ) ee = (h + ) ann , which 
allows us to solve for z as a function of the plasma temperature. We find that for the highest 
temperature of interest, T = 10 11 K, z has a maximum value ~ 0.1, and that it decreases 
monotonically in cooler plasmas. This implies that our low- z assumption is consistent. 

Given z we can compute the total number density of electrons in the disk, n e = n p (l + z), 
and number density of electrons plus positrons, n± = n p (l + 2z). 

2.3 Cooling Processes 

Since ions are considerably more massive than electrons and positrons, viscous heating will affect 
primarily the ions (SLE, Phinney 1981, Rees et al. 1982), which will then transfer this energy 
to the electrons. In their paper NY assume that the energy transfer occurs only via Coulomb 
interactions, ignoring possible non-thermal coupling mechanisms (e.g. Begelman &: Chiueh 1988). 
Since Coulomb interactions are not efficient at low densities, the ions are always hotter than the 
electrons. In this paper we assume that the ions and electrons have been brought into thermal 
equilibrium by some process, without explicitly specifying the actual coupling mechanism. 

The cooling of the hot plasma occurs mainly through the electrons and positrons. In our 
model we consider three cooling processes: bremsstrahlung, synchrotron, and Compton cooling 
by both bremsstrahlung and synchrotron photons. 

Bremsstrahlung Cooling The free-free radiation in a plasma is produced through five differ- 
ent types of interactions: electron-electron (e~e~), electron-ion (e~i), positron-electron (e~e + ), 
positron-positron (e + e + ), and positron-ion (e + i). However, e~i and e + i processes can be com- 
bined together since their corresponding cooling rates are identical, and the same is true for e~e~ 
and e + e + processes. Thus the total bremsstrahlung cooling rate per unit volume is just the sum 
of the rates for e^i, e ± e ± , and e~e + processes; and we can write % r = <Qi J r q^ e + provided 
we replace n p n e by n p n± and n 2 by (n 2 + n\) in our calculations of q~ { and q~ e respectively. 

Following NY, we adopt the expressions from Stepney & Guilbert (1983) and Svensson (1982) 
for the e ± i, e ± e ± , and e~e + bremsstrahlung cooling rates: 

q~ = 1.48 x 10 -22 n p n± F ei ergcm^s" 1 ; (7) 

_ _ f 2.56 x W- 22 (n 2 e + n 2 + )9 1 - 5 (l + 1.19 + 2 - 1.250 2 - 5 ) erg cm^s" 1 , if 9 < 1, 
qee \ 3.42 x lO- 22 (n 2 + n^)0 (In (1.1230) + 1.28) erg cm" 3 a" 1 , if > 1; [ ' 

j 3.43 x 10~ 22 n e n + (0 a5 + 1.70 2 ) ergcm^s- 1 , if < 1, 

q+ ~ I 6.84 x 10~ 22 n e n+ (In (1.1230) + 1.24) erg cm" 3 s" 1 , if > 1; ■ 



where 




96» 



(1 + 1.781 9 LM ), if 0< 1, 

(In (1.1230 + 0.48) + 1.5), if > 1. 



(10) 
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Synchrotron Cooling Synchrotron cooling is very important in our model, both because we 
assume an equipartition magnetic field in the disk and because we consider relativistic tempera- 
tures for the electrons and positrons. 

Pacholczyk (1970) has calculated the spectrum of synchrotron emission by a relativistic 
Maxwellian distribution to be 

-so ^vn± T (x M \ , _-3„-l 



where 



e s du = 4.43 x 10~ dU - * / — r dv erg cm~ 6 s" 1 , (11) 
K2{l/v) \smr 



2u eB 



XM = o — T^i u 



SuqO 2 ' 27rm e c' 

</> is the angle between the velocity vector of the electrons and the direction of the local magnetic 
field, and I(x) is a tabulated function. Positrons produce synchrotron radiation at the same rate 
as electrons do and therefore in equation (|ll]) we have used the total number density, n±, instead 
of n e . 

For an isotropic velocity distribution we can average I(xm/ sin </>) over (ft to get a new function 
I'{xm) for which a fitting function was found by Mahadevan, Narayan &; Yi (1996): 



, 4.0505 / 0.40 0.5316 \ / i/ 3 \ 
l'(x M ) = 1 + — + exp (-1.8899xif ) . (12) 



We substitute this function for I(xm/ sm.(ft) in equation (|il|). 

Equations ( |TTj ) and (|l^) are valid only for optically thin emission. However, below some 
critical frequency v c the emission becomes self-absorbed and that has to be taken into account 
in computing the total cooling rate. We estimate u c as the frequency at which the synchrotron 
emission from a thin annulus of height 2H, radius R, and thickness AR is equal to the blackbody 
emission (in the Rayleigh-Jeans limit) from the upper and lower surfaces of the annulus (NY 
presented a similar argument but they used a spherical rather than cylindrical geometry). This 
condition gives us the equation: 

v 2 

2H ■ 2-kRAR -e s dv = 2- 2-kRAR ■ 2ir-^kT du, (13) 

cr 

which we can solve numerically for given values of R, H, and T to obtain v c . 

Given the value of the critical frequency, we can estimate the total synchrotron emission per 
unit volume as follows. We assume that for frequencies below v c the emission is completely self- 
absorbed so that the volume emissivity can be approximated by the blackbody emission from 
the surface of the disk divided by the disk volume. Above v c the emission is optically thin and 
equation ( |Tl| ) can be used. To get the total cooling per unit volume by synchrotron emission, we 
then integrate over frequency: 



2ttR 2 f u c v 



2 



2-ir—kT du + I e„ dv 



V 



2HirR 2 Jo c 2 

2irkTu^ 9a n+ 

c + 6.76 x 10" 28 



3Hc 2 ' K 2 (l/0)a 1 /6 



1V2 V 2 



+ 



19/4 

a. 



(14) 
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where the parameters a±, 0,2, 03, 04 are defined as 



2 0.4 0.5316 1/3 , , 

ai = ^~02' ° 2 = ~T/4' ° 3 = T72 - ' ° 4 = L8899a l ■ I 15 ) 

Compton Cooling Since our solutions are hot, Comptonization of soft photons by hot elec- 
trons becomes an important cooling mechanism. It is an especially important process for higher 
values of the accretion rate, when the scattering depth is of order unity and synchrotron emission 
is the dominant photon production mechanism, i.e. most photons are initially very soft. 

We compute the Compton cooling rate using the Comptonized energy enhancement factor 77, 
which is defined as the average ratio of the energy of a photon at escape to its initial energy. We 
adopt the following prescription for 77, which is derived in Appendix A: 

v = e'V-V [1 - P (j m + 1, As)} + Vma X P(jm + 1, s), (16) 

where P(a,x) = [1/T(a)] Jq t a ~ 1 e~ t dt is the incomplete gamma function and 

A = 1 + 46 + W6 2 , 

S = T es + T es , 
j m = hi (rjmax)/ 111 (A), 

Vmax — • \ J- ' ) 

The parameter A is the average increase in energy of a soft photon per scattering for a Maxwellian 
distribution of electrons (and positrons) at temperature 9, T es = 2n±axH is the Thomson optical 
depth (counting both electrons and positrons), j m is the number of scatterings required for the 
maximum possible energy enhancement r] max , and v is the initial photon frequency. By definition, 
the Comptonized flux is the original flux multiplied by 77. 

Since the synchrotron spectrum is strongly peaked at the critical frequency v c , we estimate 
the cooling rate due to the Comptonization of the synchrotron radiation simply by computing r] 
for v = v c : 

Qc,s = Vfa) Q7- ( 18 ) 

The spectrum of bremsstrahlung radiation is practically flat for frequencies between the 
bremsstrahlung self-absorption frequency, z^ r , and the exponential cutoff at v = kT/h. The emis- 
sion at frequencies v < tv is negligible compared with the integrated bremsstrahlung cooling rate. 
Therefore, we can approximate free-free emission per unit frequency as the ratio — z-v)- 

We estimate as the frequency at which the free-free absorption opacity, t<JJ (which is frequency 
dependent) is equal to the electron scattering opacity, n es . Below z-v Comptonization is not im- 
portant, since absorption dominates over scattering. Therefore to calculate the Comptonization 
enhancement we perform the following numerical integration: 

kT_ -, — kT_ _ 



Qcbr= •niy) J r Ldu = *iy) kT du - ( 19 ) 

Jv br av Jv br — - v hr 

Cooling of the Optically Thick Gas The total cooling rate for the optically thin gas is 
simply the sum of the two components calculated above: 

<f = Qc,br + Qc,s- ( 2 °) 
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This rate is appropriate for the hot solutions, which are generally optically thin. However, to 
reproduce the equilibrium solution corresponding to the cool optically thick Shakura & Sunyaev 
disk, we use the generalized cooling formula given by NY: 

^ r , (21) 



l 5T + ^ + ^Tt^- M + q - 

where r = r es + T a b s is the total optical depth of the disk in the vertical direction and r a 6 s = 
(H/4<jbT 4 ) (qcbr + Qcs) 1S t ne optical depth for absorption. For a small optical depth, equation 
(PH) reduces to equation (pp|) , and therefore this formula can be used in the limit of both very 
high and very low optical depth. 

Radiation Pressure The total cooling rate allows us to compute the radiation pressure in the 
disk, 

» = V' T+ ;7s>- (22) 

The gas pressure in the accreting flow is simply: 

p = pkT (J— + J—\ , (23) 
\{J.im u id e m u J 

where \i{ and fj, e are the effective molecular weights of the ions and electrons, given by 

= TTsx = L23 ' <" = TTx = 1M ' (24) 

with numerical values computed for a hydrogen mass fraction X = 0.75 (NY). Given p r and p g , 
the parameter x, defined in §2.1, is readily computed: 

X = P/Pm=Pg/(Pr+Pg)- (25) 

3 Properties of the Solutions 

An equilibrium thermal state for the accretion flow requires that viscous heating minus the 
advected energy is exactly balanced by the total radiative cooling at each radius: 

q + (l-f) = q-- (26) 

For fixed values of m, r, m, and /3 m , the equations (p])-(p6|) comprise a closed set which can 
be solved to obtain all the relevant parameters of the accretion flow. We solve the equations 
numerically using T and x as free variables. We begin by assuming arbitrary values for T 
and x an d use equations (@)-(]6|) to calculate all the parameters of the plasma, including the 
pair fraction. Then we use equations (|7|)-(|2"4|) to calculate the cooling rate q~ and the gas and 
radiation pressures, p g &ndp r . Now we compute x = Pg/ (Pr+Pg) (see equation (p5|) ) and compare 
the result with the initially assumed value. If the two values are not equal, we vary x keeping 
T fixed until there is agreement between the initial and final x- We then vary T, optimizing x 



at each step, until the energy equation (26) is satisfied. At this point we have a self-consistent 
equilibrium solution. 
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All calculations presented below were done for an equipartition magnetic field strength in the 
disk, j3 m = 0.5. However, we have found that models with (5 m = 0.95 (where magnetic pressure 
constitutes only 5% of the total pressure) produce practically the same results. 

Figure |l| illustrates the variations with T of the two sides of equation (p6|); the solid line 
shows q + (l — f) and the dashed line shows the total radiative cooling q~ '. These results are for 
a single-temperature accretion disk of mass accretion rate M = 10~ 5 MEdd around a black hole 
of mass M = 10M Q at radial distance R = WRgchwi X has been optimized at each T. The three 
points where the two curves intersect, labeled as 1, 2, and 3 on the figure, correspond to three 
equilibrium states of the accretion disk where the energy balance equation (p6|) is satisfied. The 
equilibrium point 1 corresponds to the standard Shakura & Sunyaev (1973) thin disk; the middle 
point is the one-temperature equivalent of the unstable Shapiro, Lightman, & Eardley (1976) 
solution; the equilibrium point labeled 3 is a new hot solution which which does not seem to have 
been discussed before in the literature. It is related to the advection-dominated one-temperature 
solution of Abramowicz et al. (1995), but with some important differences as we argue later. 

Note that the slope of the cooling curve in Figure [l] changes twice. At T ~ 10 6 ' 2 K the 
change of slope indicates a transition from an optically thick to an optically thin flow. There is a 
corresponding kink in the heating curve, since at this point the cooling rate is at its maximum and 
radiation pressure briefly dominates over the gas pressure, causing a sharp decrease in gas density. 
For T < 10 6 - 2 K, the gas cools blackbody and the cooling rate increases with increasing 
temperature. For T > 10 6,2 K, the dominant cooling process is optically thin bremsstrahlung, 
which is less efficient at higher temperatures. At T ~ 10 9 ' 5 K the slope of log(g~) changes 
again, indicating a point where synchrotron radiation becomes the dominant cooling process. 
In contrast to bremsstrahlung, the synchrotron emission increases with increasing temperature, 
hence the positive slope of the cooling curve at these high temperatures. 

By inspection it is obvious that equilibria 1 and 3 are thermally stable. In both of these 
equilibria, if the gas is perturbed to a higher temperature, the rate of cooling becomes greater than 
the rate of heating, allowing the medium to cool back to the equilibrium value of T. Similarly, 
if T becomes slightly smaller, the heating rate exceeds the cooling rate and the gas heats up 
back to the equilibrium temperature. The equilibrium 2 on the other hand is unstable. Here a 
small increase or decrease in T leads to a runaway situation where T deviates progressively more 
rapidly from the equilibrium value. 

We have computed the properties of the accreting gas corresponding to the three equilibria 
for different values of M. The resulting curves for a central black hole of mass M = IOMq at 
fixed radial distance R = lORgchw are plotted in Figure |2| for a = 0.1 (heavy line) and a = 1.0 
(thin line). The curves are labeled 1, 2, and 3 to indicate correspondence with the equilibria 
labeled 1, 2, 3 in Figure |l|. The six panels in Figure || show how various physical properties of 
the solutions depend on the accretion rate. The unstable SLE solution (branch 2) is indicated 
by a dashed line. 

For a = 0.1 the topology of our results is qualitatively similar to that of the corresponding 
solutions for two-temperature disks described by Chen et al. (1995). We see that the two hotter 
branches merge together at some critical accretion rate (M cr n), above which the cooling in the 
disk is so efficient that the only equilibrium allowed is the standard thin disk solution. For 
a = 1 the hot branch extends to much higher values of the accretion rate; in fact, the value of 
Merit goes up by ~ 2 orders of magnitude, in agreement with the results of Abramowicz et al. 
(1995) and NY. Somewhat surprisingly, the topology of the solution curves for a = 1 continues 
to be similar to that for a = 0.1. This is different from the result reported by Chen et al. 
(1995), who found that for a = 1 the middle branch merged with the cold disk solution, whereas 
the hottest branch extended independently upward to become the optically-thick advection- 
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dominated solution discovered by Abramowicz et al. (1988). The value of the critical a at which 
the topology changes is, however, sensitive to the details of the radiation processes included in 
the model (M. Abramowicz, private communication). In fact, we do find the same topology 
change in our results when we push a above ~ 1.7. 

For both a = 0.1 and a = 1, the curves corresponding to the three equilibria in Figure §(a) 
have positive slopes, i.e. dM /dT, > 0. This implies that all three solutions are viscously stable 
(cf. Frank, King, & Raine 1992, p. 103). Note also that our new hot solution is optically thin 
for all values of M where the solution exists (Figure ||[b] ) and cools mainly through synchrotron 
and Comptonized synchrotron processes (Figure ^[c] ) . 

Figure ^ directly compares the two-temperature (thin line) and one-temperature (heavy line) 
equilibrium solutions for a disk with a = 0.1. For temperatures below T ~ 10 9 K the two 
curves merge, since the electron and ion temperatures in the two-temperature solution are nearly 
equal (Figure ^[c]) and therefore, two-temperature and one-temperature solutions are effectively 
the same. However, for T ^ 10 9 K the two solutions differ significantly. The most important 
distinction is in the value of the advected energy fraction / (Figure |3][b]). For the hot two- 
temperature solution, we have / — ► 1 at all M, showing that this branch is always advection- 
dominated. In contrast, we see that our one-temperature hot solution is cooling-dominated 
(/ < 0.1) for a wide range of accretion rates. This is a completely new result, as prior to this 
the only hot stable solutions known were all advection-dominated (NY, Abramowicz et al. 1995, 
Chen 1995, Chen et al. 1995). Since the cooling rate is given by q~ = q + (l — /), we see that, at 
a given accretion rate, our single-temperature disk with its low value of / is considerably more 
luminous than the equivalent two-temperature flow. 

In the case of the two-temperature advection-dominated hot solutions discussed by NY and 
Abramowicz et al. (1995), the thermal stability of the solution to long- wavelength perturbations 
has been shown to be the direct result of advection. However, our hot one-temperature solution 
branch 3 is cooling-dominated and has very little advection. Why is this solution stable? The 
answer is obvious from Figure |l], where we see that the introduction of synchrotron cooling causes 
the cooling curve q~(T) to be steeply positive for T > 10 9 5 K. Therefore, thermal stability is 
achieved in this solution because of the usual reason, namely that the cooling increases more 
rapidly than the heating. Note that we would not have obtained this result if we had included only 
bremsstrahlung cooling. In fact, Abramowicz et al. (1995) considered the pure bremsstrahlung 
one-temperature case and they obtained a hot advection-dominated solution. Thus, the inclusion 
of synchrotron cooling has an important effect on the basic physics of the accretion flow. Although 
a full linear stability analysis along the lines of Kato et al. (1995) has not been done on our 
new solutions, we suspect that these solutions are stable to all linear perturbations, regardless of 
wavelength. 

Figure |] illustrates the radial structure of a disk with a = 0.1 for a fixed value of the accretion 
rate, M = W- 4 M Edd . The inner edge was chosen to be at R = 3Rg c hw: corresponding to the 
last stable orbit. As expected, we see that the two-temperature (thin line) and one-temperature 
(heavy line) solutions merge at large radii where the temperature is low and even the two- 
temperature plasma has effectively a single temperature (see Figure §[c]). Thus, away from the 
center, the flow is advection-dominated (see Figure ||[b]) and almost spherical with H/R ~ 1 
(Figure ||[d]), regardless of whether we have a two-temperature or one-temperature solution. 
However, closer to the accreting object, the flow becomes increasingly cooling-dominated in the 
one-temperature case and settles into a disk with H/R ~ 0.15 at the inner edge. The resulting 
geometry - thinnish disk near the inner edge and quasi-spherical flow at larger radii - is unusual 
for accretion disk models. Another interesting feature of our one-temperature solution is that the 
temperature profile of the disk is not monotonic but has a well defined maximum at R ~ AORschw- 
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Therefore, the hottest part of the spectrum does not come from the inner edge of the disk. 

NY showed that their advection-dominated hot two-temperature solution exists only for ac- 
cretion rates M less than a critical rate Merit i which depends on a and R (see also Abramowicz et 
al. 1995). We find a similar result also in the case of our cooling-dominated hot one-temperature 
solution. Figure [B| shows the critical accretion rate M cr n for our solution as a function of R as 
well as the variation of the gas temperature. The two panels on the left (Figure 5[a,b]) compare 
the results for the one-temperature and two-temperature solutions for a = 0.1. We see that the 
electron temperatures for the two solutions are nearly equal, but the behavior of M CT n is quite 
different in the two cases. For a two-temperature flow, M cr a is roughly independent of R for 
R 5s 10 3 Rschw and is given approximately by M ~ 0.3a 2 Mem (see NY and Abramowicz et al. 
1995 for a discussion of the a 2 scaling). In contrast, M cr n decreases almost linearly with de- 
creasing R in the one-temperature case, so that M ~ 310 _3 a 2 (i?/ RscHw)Mem for R 10 3 Rschw 
Thus, near the black hole, M cr n for the one-temperature solution is smaller than that of the 
two-temperature solution by as much as two orders of magnitude. 

The steep decrease of M cr u with decreasing R leads to an interesting effect which is already 
seen in Figure |j. For a wide range of M, it is possible to have M < M cr u at large radii, but to have 
M > M cr n at small radii, with a transition at some critical radius R m in- For the example shown 
in Figure |], we have R m in = SRschw in this case, we can have a hot one-temperature solution for 
R > R m in, but the solution disappears at R = R m i n - Below R m in, the only equilibrium solution 
available is the cool thin accretion disk. (Note that, because of the very low M ~ W~ 4 MEdd, 
the Lightman-Eardley instability does not affect the thin disk solution in this case.) Thus, the 
example shown in Figure || has an extremely unusual structure. For 3Rs c hw < R < SRschw, there 
is a standard thin disk. Then, for 5Rs c hw < R < ^QRschw, we have a hot, but cooling-dominated 
thickish disk. Finally, for R > 40Rs c hw, we have a regular advection-dominated quasi-spherical 
flow. There are several discussions in the literature of models where an outer thin disk makes 
a transition to a hot inner flow (e.g. SLE, Wandel & Liang 1991, Narayan, McClintock, & Yi 
1996, Narayan 1996), but here we have an example of a model which does the opposite, namely 
transforms from an outer hot advection-dominated flow via a cooling-dominated hot flow to an 
inner thin disk. (See Melia 1994 for the closest example discussed previously in the literature.) 

Figure |(| shows how the properties of the flow change as the accretion rate is reduced. There 
are some interesting effects. First, with decreasing M, the critical radius R m in disappears so that 
there is no longer any reason for the gas to make a transition to a cool disk on the inside. As M 
decreases still farther, even the cooling-dominated zone of the hot disk disappears, and advection 
becomes more important throughout the flow. Consequently the disk thickens considerably, and 
the maximum in the temperature profile at large radius becomes less apparent. For a = 0.1 and 
M / Mfidd < 10~ 8 the flow is advection-dominated and nearly virial at all radii. 

Figure |6|c] shows the relative importance of Comptonized synchrotron cooling as a function 
of R at various M. We see that synchrotron and Comptonized synchrotron cooling is dominant 
in the inner regions, R < lOORschw, regardless of the accretion rate. In the outer disk, on the 
other hand, where the flow is cooler and more advection-dominated, the cooling is mostly due to 
bremsstrahlung emission. 

Though one-temperature solutions are restricted to smaller values of M than the correspond- 
ing two-temperature solutions, nevertheless we should note that for any given M they are signif- 
icantly more luminous. This is because these flows are cooling-dominated at small radii where 
most of the radiation is emitted. Figure 0(a) illustrates this luminosity increase in the case of 
a = 0.3. Here we plot the total disk luminosity as a function of M for one-temperature (solid line) 
and two-temperature (dashed line) solutions. Note that at M = Mem the one-temperature 
disk is ~ 10 times more luminous than its two-temperature counterpart. 



11 



Figure 0(b) shows how the total luminosity of our one-temperature solutions changes with 
a. For a given value of M, higher-a solutions are less luminous, since they are more advection- 
dominated (see Figure @[e]). However, higher values of a allow the hot solution to survive up 
to higher accretion rates (see Figure |5j[c]) with correspondingly higher luminosities. Thus, the 
most luminous system that we can model as a hot one-temperature disk extending down to 
R = 3Rschw requires a = 1.0 and M = W~ 3 MEdd, and has a total luminosity L = W~ 28 LEdd- 
We could in principle use higher values of M, but then as discussed above, we will have to include 
a cool thin disk at small radii (R < R m i n ). 

4 Discussion 

Most of the work done so far on modeling hot accretion disks has been based on two-temperature 
flows, where the electrons are much cooler than the ions (SLE; Rees et al. 1982; Melia & Misra 
1993; Narayan &; Yi 1995b; Narayan et al. 1995, 1996). However, if there exist mechanisms 
of energy transfer between ions and electrons that are more efficient than Coulomb scattering 
(Phinney 1981, Rees et al. 1982, Begelman &: Chiueh 1988), the plasma might be well-coupled 
and the ions and electrons may have the same temperature. In this paper we have investigated 
how the properties of the hot, stable two-temperature solution described by NY and Chen et al. 
(1995) are altered if ions and electrons in the accreting plasma are constrained to have the same 
temperature. 

We start from the set of equations developed by NY, adding radiation pressure and particle- 
particle pair production processes, and improving some of the radiative cooling calculations. We 
solve these equations for various choices of the viscosity parameter, a, and magnetic pressure 
parameter, /3 m , explicitly setting the ion and electron temperatures equal to each other. The 
resulting three branches of solutions which we obtain (Figure [|) correspond to the standard 
thin disk (Shakura & Sunyaev 1973, Novikov &: Thorne 1973, Lynden-Bell & Pringle 1974), 
the one-temperature equivalent of the unstable hot solution discovered by SLE, and finally a 
new thermally stable hot solution whose properties we discuss in detail (§3). The new solution, 
which is identified by label 3 in Figures 1 and 2, is the one-temperature equivalent of the two- 
temperature advection-dominated solutions discussed by NY (see also Abramowicz et al. 1995, 
Chen et al. 1995). 

For a viscosity parameter a = 0.1, we find that the topological relationships among the 
three solution branches are similar to those found by Chen et al. (1995) and NY for the two- 
temperature solutions. Specifically, we see that the two hot branches merge at some maximum 
value of the mass accretion rate, which we call M cr a, such that for M > M cr u the only solution 
available to the flow is the standard thin Shakura & Sunyaev disk configuration. At higher a, we 
confirm the result of Chen et al. (1995) that the topology changes. However, the critical a where 
the change occurs is 1.7 (at R = WRschw) for the one-temperature case we consider, rather then 
0.3 in Chen et al. (1995). 

We have computed the value of the critical mass accretion rate M cr u as a function of radius for 
different values of a (Figure |5|[c]) and P m . We find that M C rit varies roughly as a 2 (Abramowicz 
et al. 1995, NY) and is relatively independent of P m . A new result is that M cr n varies almost 
linearly with R for R <> Rschw, since in a two-temperature flow M cr u is roughly independent of 
R in that region. The maximum value of the accretion rate at which the hot solution exists near 
the inner disk radius (R = 3Rs c hw) is M cr i t ~ 3 x W~ 3 a 2 MEdd- In comparison, two-temperature 
hot solutions have M C rit ~ 0.3a 2 MEdd- 

Though our single-temperature solutions are limited to lower mass accretion rates than those 
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allowed by two-temperature models, this is compensated to some degree by strong cooling in the 
inner regions, where / < 0.1 for M > 10 -3 M cr i t (3Rs c hw) (Figure ||[f]). As a result, over the 
range of mass accretion rates, 10 -3 < M / M cr it(3Rs c hw) < 1 5 our solutions are significantly more 
luminous than the corresponding two-temperature solutions (Figure 7[b]). 

This work represents the first example of a hot thermally stable accretion solution which is 
also cooling-dominated and therefore an efficient radiator. Because advection plays a very minor 
role in our new solutions, the thermal stability is not the result of advection as in the equivalent 
two-temperature solutions (NY, Abramowicz et al. 1995) or in the pure bremsstrahlung one- 
temperature solutions (Abramowicz et al. 1995). We show that the thermal stability of the hot 
one-temperature solution is primarily the result of including synchrotron emission, which leads 
to a rapid increase of the cooling rate with increasing temperature (see Figure [l]) . 

An interesting feature of our models is the unusual geometry of the flow. Since M cr n decreases 
with decreasing R in the inner parts of the disk (see Figure ||) , for a wide range of M the disk has 
a very unusual structure. In the inner regions, where M > M cr it(R), the only possible solution 
for the accreting gas is the Shakura & Sunyaev thin disk. Note that since we are restricted 
to relatively low mass accretion rates, the thin disk is viscously stable even at R ~ 3Rs c hw 
Further out, the hot solution exists, but it is cooling-dominated and therefore has a disk-like 
geometry. Finally, in the outer parts of the disk the flow is fully advection-dominated and 
becomes nearly quasi-spherical. This interesting result, however, presents a problem for our 
calculations. The inner radiative zone of the disk is cooler than the outer advection-dominated 
zone, so the relatively soft photons emitted on the inside will Compton cool the hot outer gas in 
a way that is incompatible with the local treatment of the radiative processes which we have used 
in this paper. The effects of the non-local radiation field need to be incorporated self-consistently 
into future calculations. 

Since the gas temperatures in our solutions are on the order of T ~ 10 10 K, pair production 
and annihilation processes must be taken into account. In our calculations we have included 
particle-particle pair production rates (White & Lightman 1989) which give values for the pair 
fraction z less than a few percent, but we neglected photon-photon and photon-particle pro- 
cesses. We believe that this treatment is justified for low values of m (low optical depth), since 
bremsstrahlung radiation and Comptonization become unimportant and the cooling is dominated 
by soft synchrotron photons which are incapable of producing pairs. However, for M ^ 0.1M cr n, 
Comptonization of synchrotron radiation begins to dominate and photon processes may become 
important. To place a limit on these processes we computed the equilibrium pair density tak- 
ing into account photon-photon and photon-particle pair production in the Wien peak of the 
Comptonized spectrum, using the rates derived by Svensson (1984). The resulting values of z 
are significantly higher than the values computed taking into account only particle-particle pair 
production, as M approaches M cr n. However, even at M = M cr u, the maximum value of z 
that we find with the inclusion of photon interactions is only ~ few%. Kusunose (1996) and 
Abramowicz (1996) confirm that pair processes are unimportant in the hot optically thin flows 
which they have studied. Thus, we believe that our treatment of pair processes is reliable in the 
entire parameter space we have considered. 

This work was supported in part by NSF grant AST 9423209 to the Center for Astrophysics. 
RN thanks the Institute for Theoretical Physics (NSF grant PHY 9407194) for hospitality. IY 
acknowledges financial support from SUAM Foundation. 
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A Comptonization 



Consider a point at an optical depth r es inside the medium. The mean number of scatterings of 
a soft photon escaping from this point is (Rybicki & Lightman 1979) 



r 2 

' es' 



(27) 



Let the average photon energy change per scattering be A, i.e. the initial and final photon 
energies, Ei n and Efi n , are related on average by Efi n = AEi n . For a thermal distribution of 
electrons with temperature 9 = kT/^rrieC 2 ), we have (Rybicki & Lightman 1979) 



A m 1 + 49 + 169' 



(28) 



The Compton y-parameter is simply y = s (A — 1) = s (49 + 169 2 ). 

After j scatterings, the energy gain of the photon is A J so long as Ej in is not saturated to 
the Wien regime (Efi n ~ 3kT). The maximum j corresponding to saturation is thus 



hitman/ hi A, where rj r , 



VWi 



3kT 



(29) 



The energy enhancement factor is then simply the average energy gain of the photon: 



v = Y, Ajp i+ Aim E p i> 

j=0 J=in»+1 



(30) 



where Pj is the probability that a photon will suffer exactly j scatterings. 

To compute 77, we need an expression for Pj. The probability that a photon travels a distance 
characterized by an optical depth r es without scattering is e~ Tes . To compute Pj , the assumption 
is usually made that in a medium with r es 1, successive scattering events can be treated 
independently (e.g. Dermer, Liang, & Canfield 1991), so that Pj can be written as e~ Tea (l — 
e~ Tes y . However, for spherical geometry the scattering probability decreases after successive 
scatterings, because the mean square distance of a photon from the center increases with each 
scattering. Allowing for this accurately requires Monte Carlo calculations. One limit however, 
is straightforward, namely the case when the photon continues to travel radially outward after 
each scattering. In this case, the probability Pj is given by the Poisson formula 



(31) 



With the above expression for Pj we can write rj as 



^ (Asy . , 

.. + A 3r ' 

3=0 



i! 



j=jm + l 



E 



(Asy 



+ A ' m E S ~r 



To evaluate the sums in equation (p 



we define the following function: 

00 j 

m= E fr- 

j=jm+l J ' 



(32) 



(33) 
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Using integration by parts and the fact that /(0) = 0, it is easy to show that 



/OK 



1 



e y y- ,ra dy 



(34) 



and therefore, 

/(*)= E ^ = —y\ e-y y ^dy = e x P(j m + l,x), (35) 

where P(j m + 1, aj) is the incomplete gamma function. Substituting this result into equation ( |32| ) 
yields the final result: 
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_ e ^p(j m + l, As)) + A^e s P(j m + 1, s) 



VrnaxP(jrn ~i~ 1) s). 



(36) 
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Figure 1: The rates of heating (solid line) and cooling (dashed line) of the accreting single- 
temperature gas plotted vs. gas temperature, T. Intersection points of the two curves mark 
equilibrium states of the system. Point 1 corresponds to the standard thin disk solution, point 2 
is equivalent to the unstable SLE solution, and point 3 is our new, hot, thermally and viscously 
stable, solution. 
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Figure 2: Thermal equilibria of single-temperature accretion disks with (3 m = 0.5 around a 
WMq black hole at a radius of R = WRschw The heavy and thin lines represent solutions with 
a = 1 and a = 0.1 respectively. The three branches labeled as 1, 2, and 3 correspond respectively 
to the standard thin disk, unstable SLE disk (indicated by a dashed line), and our hot and stable 
solution (as in Figure 1). In the six panels we plot (a) surface density, £ (gem -2 ), (b) total 
optical depth, (r es + T a b s ), (c) fraction of the total cooling due to Comptonized synchrotron 
radiation, (d) temperature, T (K), (e) the vertical scale height, H/R, and (f) advected energy 
fraction, /, all as functions of the accretion rate M. Note that the hot solution branch (labeled 3) 
is optically thin, has essentially spherical geometry (H/R « 1), and cools mainly by synchrotron 
radiation. 
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Figure 3: Thermal equilibria of single-temperature (heavy line) and two-temperature (thin line) 
accretion disks with (3 m = 0.5 and a = 0.1 around a 10M Q black hole at a radius of r = 10. 
The thin disk branch is not shown. The four panels show (a) surface density, £ (gem -2 ), (b) 
advection parameter, /, (c) temperature, T (K), and (d) the scale height, H/R, plotted as 
functions of M. The unstable branch is shown by dashed lines. On the stable branch, the single- 
temperature solution is less advection-dominated and has a higher electron temperature than the 
two-temperature solution. 
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Figure 4: Radial profiles of thermal equilibria for a single-temperature accretion disk with 
An = 0.5, a = 0.1, and M / MecLcL = 10~ 4 around a 10M Q black hole (heavy line). For comparison 
we also plot the hot stable branch of the two-temperature disk (thin line). As expected, the 
two models give identical results at large radii. Closer to the black hole, however, the single- 
temperature solution becomes cooling-dominated while the two-temperature solution remains 
advection-dominated. Note the maximum in the disk temperature at R « 40Rs c hw 
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Figure 5: Critical accretion rate M cr n (upper panels) and corresponding temperature (lower 
panels) are plotted as functions of the radius. Panels (a) and (b) compare the results for the one- 
and two-temperature models for (3 m = 0.5 and a = 0.1. The critical M differs significantly in the 
two models, but the gas temperature in the single-temperature model is practically always equal 
to T e in the two-temperature model. Panels (c) and (d) show the results for a single-temperature 
disk with (3 m = 0.5 and two different values of a. Note that although a = 1.0 allows models 
with values of M C rit higher by two orders of magnitude, the corresponding temperatures do not 
depend strongly on a. 
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Figure 6: Radial profiles of hot single-temperature solutions, labeled by the values of the 
accretion rate in units of Me<u- The other parameters have the same values as in Figure 4. Note 
that for M = W~ A MEdd the hot branch does not extend all the way to the inner edge of the 
disk at R = SRs c hw Therefore, the only stable configuration at very small radii is the thin disk 
solution which is also shown on the figure. 
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Figure 7: (a) Compares the integrated luminosities for one- and two-temperature hot accre- 
tion disks with (3 m = 0.5, a = 0.3, and M = 1OM . The single-temperature disk is clearly 
more efficient than a two-temperature disk for a given M since its interior regions are cooling- 
dominated. However, the limiting M cr u of the one-temperature solution is lower than that of 
the two-temperature flow by two orders of magnitude, (b) Integrated luminosity of a hot single- 
temperature accretion disk around a 10M Q black hole as a function of the mass accretion rate, 
for j3 m = 0.5 and three different values of a. Each curve was calculated up to the maximum 
value of the accretion rate given by M max = M cr i t {R = 3Rs c hw)- The dotted line shows the 
dependence expected if the radiative efficiency is 10%. Note that disks with lower values of a 
are more luminous at a given M, but have a lower limiting accretion rate Merit- 
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